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We study the ground state properties of a quantum antiferromagnet on the kagome lattice in 
the presence of a magnetic field, paying particular attention to the stability of the plateau at 
magnetization 1/3 of saturation and the nature of its ground state. We discuss fluctuations around 
classical ground states and argue that quantum and classical calculations at the harmonic level do 
not lead to the same result in contrast to the zero-field case. For spin S — 1/2 we find a magnetic 
gap below which an exponential number of non-magnetic excitations are present. Moreover, such 
non-magnetic excitations also have a (much smaller) gap above the three-fold degenerate ground 
state. We provide evidence that the ground state has long-range order of valence-bond crystal type 
with nine spins in the unit cell. 

PACS numbers: 75.10.Jm, 75.60.Ej, 75.45.-|-j 



I. INTRODUCTION 

The appearance of exotic quantum phases in systems 
described by two-dimensional frustrated antiferromag- 
nets is presently the subject of intense research (see e.g. 
Refs. 1,2 for recent reviews). The Heisenberg antifer- 
romagnet on highly frustrated lattices such as the py- 
rochlore and kagome lattice has a huge degeneracy of the 
classical ground state such that no magnetic order arises 
at any temperature (see e.g. Ref. 0for a recent review). 
At the quantum level one may then obtain different exotic 
phases without magnetic (Neel) order. One such phase 
is the so-called 'valence-bond crystal' which is character- 
ized by formation of local singlets in a long-range ordered 
pattern. An even more exotic phase, namely one without 
any kind of long-range order, is suspected to arise in the 
S = 1/2 Heisenberg model on the kagome latticei^i^*Sii 
In the latter case, there is a small spin gap and, although 
this is still under discussion^ the ground state is sus- 
pected to be disordered. In particular, a huge number 
of singlets (exponentially growing with the system size) 
is found inside the spin gap which are reminiscent of the 
classical degeneracy. 

The spin 5=1/2 kagome Heisenberg antiferromag- 
net (KHAFM) is realized e.g. in volborthite,^° although 
presumably in some distorted form. Another possible 
realization is given by atomic quantum gases in optical 
latticesiii In the latter case, magnetization corresponds 
to particle number and a magnetic field to chemical po- 
tential, opening the possibility to perform experiments 
for the behavior of the spin model in a magnetic field. 

The magnetization process of the KHAFM has been 
studied theoretically both for classicaliSiiii^ as well as 




FIG. 1: Boundaries of the (M) = 1/3 plateau as a function 
of the anisotropy A for difi'erent lattice sizes (see legend) and 
the transition to saturation (M) = 1 for the thermodynamic 
limit (thin full line). Inset: kagome lattice with an ordered 
state of the valence-bond crystal type at (M) = 1/3: circles in 
certain hexagons indicate local resonances between different 
Neel configurations on the hexagons, arrows indicate spins 
which are aligned with the field. 

quantum spins ni^^i^iiSiii Numerical results for the mag- 
netization curve of the 5=1/2 Heisenberg model exhibit 
among others a clear plateau at 1/3 of the saturation 
magnetizationi^iiSiiSiii (see also Fig.^. For the classical 
KHAFM at one third of the saturation field thermal fluc- 
tuations select coUinear states, but there appears to be 
no real order li^ For the 5=1/2 KHAFM we will argue in 
this paper that the state with magnetization (M) = 1/3 
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FIG. 2: Low-lying excitations above the (M) = 1/3 plateau 
for the S = 1/2 Heisenberg antiferromagnet on the A'^ = 36 
kagome lattice, (a) Full lines show all excitations with = 6 
in the given energy range, bold dashed lines the lowest excita- 
tions with = 5 and 5*^ = 7 as a function of magnetic field 
h. (b) Excitation energy versus number of states with 5*^ = 6 
below that energy. One observes a total of 100 states below 
the magnetic gap in the middle of the (M) = 1/3 plateau 
(corresponding to the largest gap to magnetic excitations). 

exhibits order of the valence-bond crystal type (the spin- 
spin correlation functions are short-rangedi^) although it 
shares some similarities with the case (M) — 0. 

In the present paper we study the XXZ model in a 
magnetic field h 

where indicates nearest neighbors on the kagome 

lattice (see inset of Fig. QJ, sf are spin-half operators 
acting at site i and A is the XXZ-anisotropy. 



II. EXACT DIAGONALIZATION FOR THE 
5" = 1/2 HEISENBERG MODEL 

First we present numerical results which have been ob- 
tained by Lanczos diagonalization of the Hamiltonian 
using the program package spinpack^ All numerical 
computations have been performed on lattices with N 
sites subject to periodic boundary conditions. 

The main panel of Fig. shows the boundaries of the 
fully polarized state (which we normalize to (M) = 1) 
and a state with (M) = 1/3 in the XXZ model (0. Fig. 
|2Ia) shows the low-lying excitations above the (M) = 
1/3 ground state computed by exact diagonalization for 
= 36 and A = 1. The dashed lines show the gap to 
states with S*^ = 5 and 7 which vanishes at the bound- 
aries of the plateau. The maximum magnetic gap occurs 
in the middle of the plateau where these two lines inter- 
sect and we will use this as a definition of the magnetic 
gap. Horizontal straight lines denote states with S*^ = 6 



[i.e. (M) = 1/3) and correspond to non-magnetic exci- 
tations. The large number of non-magnetic excitations 
below the magnetic gap is reminiscent of the classical 
degeneracy. The shape of the integrated density of non- 
magnetic excitations (see Fig.|2b)) is very similar to the 
corresponding integrated density of singlets at {M) — 
(see Fig. 2 of Ref. ^. In combination with the disor- 
dered classical ground statCj^^ one might be tempted to 
take this as evidence that also the ground state of the 
S =1/2 KHAFM at (M) = 1/3 is disordered. However, 
we will argue next that here classical and quantum fluc- 
tuations are in fact not equivalent at the harmonic level 
and then present evidence in favor of an ordered state for 
S = 1/2. 

HI. FLUCTUATIONS AROUND THE 
CLASSICAL GROUND STATE 

Classical (thermal) fluctuations were studied in Refs. 
such that we make only a few comments valid for 
non-zero magnetization and arbitrary anisotropy A. As 
was shown explicitly for A = 1 in Ref. 0i thermal fluc- 
tuations select coUinear 'up-up-down' (UUD) conflgura- 
tions at (A/) = 1/3 against the other non-coUinear con- 
flgurations that also minimize the classical energy, but 
all UUD conflgurations have the same spectra of har- 
monic fluctuations. Indeed, a local change of variables 
shows that the covering-dependent Hamiltonians of clas- 
sical Gaussian fluctuationsii are equivalent. 

The role of quantum fluctuations is however radically 
different. Now quantum commutation relations have to 
be preserved and the change of variables used for the 
classical case is no longer possible. To analyze this in 
more detail, we compute the zero-point contribution to 
the ground state energy at (M) = 1/3 for two different 
coverings with a g = and a a/S x \/3 structure, respec- 
tively (the latter can be found e.g. in Fig. 1 of Ref. 12). 
By writing the spin operators on each site in terms of 
bosonic creation and annihilation operators: 

(2) 

we obtain the Hamiltonian: 

H = Ha + ^{H2 + 0{l/^/S)), (3) 

where H2 is quadratic in creation and annihilation oper- 
ators and the 0{\/y/S) part contains higher orders. By 
Fourier transforming, we obtain: 

k ^ ^ V \ k / 

where are 3x3 and 9x9 matrices for the q = and 
the -\/3 X \/3 states since these coverings have 3 and 9 
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sublattices, respectively. No further change of variables 
is possible here since the commutation relations of the 
su{2) algebra of the spins have to be preserved. At A = 1 
one finds for the zero-point fluctuations ^ J2k^k ~ J^l"^ 
and « 0.5643 JS for the q = and the -s/S x VS state, 
respectively, demonstrating the inequivalence of the dif- 
ferent coverings at the quantum level. 



IV. EFFECTIVE MODEL FOR THE ISING 
LIMIT 

Let us now return to the extreme quantum case S = 
1/2 and study the anisotropic XXZ limit. For A ^ 1, 
we can adapt the analysis of Refs. l2lll23 of the Ising model 
in a transverse field to the XXZ model by replacing the 
expansion in the transverse field with an expansion in 
powers of 1/A. 

In the Ising hmit A = oo, the ground states are those 
states where around each triangle two spins point up and 
one down. This ground-state space of the Ising model can 
then be taken as configuration space for a perturbative 
treatment of the XF-part of the XXZ Hamiltonian, 

For sufficiently large lattices, the lowest non-trivial or- 
der is third order, flipping simultaneously pairwise an- 
tiparallel spins around a hexagon. This is described by 
an effective HamiltonianSLSS 



2cff. = A 

hexagon i 



(5) 



where the spin operator s"^ operates at the jth site 
around hexagon i and A = 3J/(2A^). Note that in 
the Ising-basis the effective Hamiltonian has only off- 
diagonal matrix elements of size A. 

The configurations of the Ising model can be mapped 
to dimer coverings of the dual lattice which in the case 
of the kagome lattice is the hexagonal lattice. Now one 
can use known results for dimer coverings22i24 to write 
down the asymptotic growth law for the number of Ising 
configurations A/'conf. on an iV-site kagome lattice: 



Konf. oc (1.11372781...) 



N 



(6) 



Exploiting results for the related quantum dimer 
model on the hexagonal lattice^^^ Moessner and Sondhi 
concluded^^ that the ground state of the effective Hamil- 
tonian is of the valence-bond crystal type. To be 
more precise, the case studied in Refs. >22tt25i corresponds 
to A < whereas we have A > 0, but there exist unitary 
transformations which change the sign of A.-^^ Hence the 
spectra of the effective Hamiltonian ^ are invariant un- 
der A —A. The three-fold degenerate ground-state 
wave functions are sketched in the inset of Fig. ^ Circles 
in one third of the hexagons denote resonances between 
the two different Neel states on the surrounding hexagon; 
a background of the remaining third of all spins points 
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FIG. 3: Main panel: Spectra of the effective Hamiltonian 
for A ^ oo with N = 36, 54, 81, 108 and 144. Inset: Scaling 
of the energy of the Ith excited state with inverse system size 
1/A^ for some selected levels. 



in the direction of the field. Note that these wave func- 
tions were argued in^^-^^ to yield a qualitatively correct 
description, but they should not be used for a quantita- 
tive analysis. Furthermore, we emphasize that due to the 
resonances, these wave functions are of a purely quantum 
nature and have no counterparts as unique states of the 
classical Heisenberg model. 

According to the above, at large A the (M) = 1/3 
state of the XXZ model on the kagome lattice should be 
three-iold degenerate with a gap to the next non-magnetic 
excitations. To check this conclusion and compare it to 
Fig. 12 let us look at the spectrum of the effective Hamil- 
tonian This effective model has a substantially re- 
duced Hilbert space (e.g. for A^ = 36 there are only 120 
states). We can therefore go to larger lattice sizes than 
in the full model. Results for kagome lattices with up to 
N — 144 sites are shown in Fig. O Additional short cy- 
cles wrap around the boundaries of the lattice for A^ < 27 
and lead to non-generic ground states of Hcs.. Accord- 
ingly, systems with A^ < 36 should not be considered and 
are not included in Fig. |31 

Two features are apparent in Fig. Oat least for the two 
biggest system sizes (A^ — 108 and 144). Firstly, there 
are two further levels above the ground state. The finite- 
size dependence of the second excited state is shown by 
I — 2 in the inset of Fig. O and indicates that it is con- 
verging to SE — > which is consistent with the expected 
three-fold degeneracy of the ground state in the thermo- 
dynamic limit. Secondly, there is a huge density of states 
emerging for SE > 1.3A. The finite-size behavior of the 
I — 5 level in the inset of Fig. O indicates that a gap of 
the order ~ 1.2 A to these higher excited states persists 
in the thermodynamic limit. For A^ ~ 108 and 144 there 
are two further levels in between. It is difficult to extrap- 
olate their energies to A ^ oo, but the behavior of the 
I = 3 excited level in the inset of Fig. O at the largest 
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FIG. 4: Overlap between the ground state of the full XXZ 
model and the effective Hamiltonian with A'^ = 36 for different 
values of the anisotropy parameter A (full line) . For compari- 
son we include the corresponding result for a different model, 
namely the = 36 triangular lattice— (dotted line) which 
exhibits a sharp drop around A = 0.77. 



values of N suggests that they retain a finite gap in the 
thermodynamic limit. Inspection of the wave- functions 
indicates that these additional low-lying levels may arise 
from the three classical a/3 x \/3 configurations. 



V. FROM THE ISING LIMIT TO THE 
HEISENBERG MODEL 

Although the effective Hamiltonian leads to higher de- 
generacies of some excited states, the TV = 36 curves in 
Figs. 13 and O have a very similar shape which can be 
taken as a first indication that the same scenario as for 
A ^ 1 also applies to A = 1 . Comparison of the overall 
scales leads to an estimate for the gap in the (M) = 1/3 
sector at A = 1 of about 0.04 J. Furthermore, the to- 
tal number of Ising configurations is very close to the 
number of non-magnetic excitations below the magnetic 
gap for A = 1 at a given system size (see also Ref. IT^ . 
Hence, the growth law © yields a good approximation 
also to the number of non-magnetic excitations in the 
Heisenberg model (A = 1). 

It is instructive to compute the overlap of the wave 
function of the full XXZ model, |full XXZ), with the 
ground state wave function, jeffective), of the effective 
Hamiltonian with the same number of spins N. The 
analysis of the effective Hamiltonian implies that one 
should study only sizes which are multiples of 9 and that 
TV = 36 is the smallest which is representative of the 
general case. However, = 36 is the biggest system 
where we have been able to study the full XXZ model 
and hence the only case we can discuss. Results for the 
overlap |(effective|full XXZ)\ are shown by the full line 
in Fig. ^ We observe that this overlap tends to 1 for 
large values of A, as expected. Furthermore, the overlap 



remains appreciable even close to the Heisenberg model 
(|(effective|full « 0.22 for A = 1), in particu- 

lar if one considers that the dimension of the symmetry 
subspace under consideration is of the order 10^. Note 
further that an analogous computation of the overlap of 
the ground state of the Ising and the full XXZ model at 
(M) — 1/3 on the triangular lattice leads to a sharp drop 
at A w 0.76 for a fixed N (see dotted curve in Fig. 0] for 
N — 36), signaling an instability of the plateau state>ii 
No such sharp drop is observed on the kagome lattice 
(full line in Fig. which we take as a sign of absence of 
phase transitions between A = cxd and « in the XXZ 
model on the kagome lattice at (M) — 1/3. In particular, 
A = 1 and oo should belong to the same phase. 



VI. DISCUSSION AND CONCLUSIONS 

To conclude, we have analyzed the low-energy spec- 
trum of the kagome XXZ S =1/2 model at magnetiza- 
tion (M) = 1/3. While the existence of a magnetization 
plateau is clear, the nature of the non-magnetic excita- 
tions over the ground state is more difficult to clarify. 
We have argued by different techniques that the ground 
state has an order of the valence-bond crystal type, i.e. 
the ground state is three-fold degenerate and there is a 
small gap to all higher excitations. While in the case 
A ^ 1 this scenario is derived from a mapping to an ef- 
fective HamiltoniaufSLSSiS^ our numerical data indicates 
that it persists down to the isotropic limit A = 1. 

One of the key differences between the present case and 
(M) = lies in the unrenormalized classical thermal and 
quantum fluctuations. In the absence of a magnetic fleld, 
they are equivalent regarding the lifting of degeneracy of 
configurations with soft modes (planar configurations). 
However, for (M) = 1/3, classical thermal fluctuations 
select the collinear UUD conflgurations and the weight 
in the free energy of any UUD covering is equivalent at 
the harmonic level. Because of commutation relations 
that have to be preserved at the quantum level, the zero- 
point corrections over the UUD configurations are not 
any more equal. This is the first indication that a spin 
liquid phase is less likely to appear than for the (Af) = 
case. 

For S* = 1/2 and A = 1 we find, for (M) = 1/3, 
an exponential number of non-magnetic excitations be- 
low the magnetic gap which are reminiscent of the clas- 
sical degeneracy; just as for {M) = Oi^ In the latter 
case the macroscopic number of non-magnetic excitations 
has been taken as evidence for a completely disordered 
ground state. Here, however, we find evidence for a fur- 
ther small gap, separating the continuum of states from a 
ground state, which has long-range order of valence-bond 
crystal type. We remark that the A^ = 36 spectrum^ sug- 
gests that candidates for valence-bond ordered states for 
(M) = would have a larger unit cell than the state 
above. Hence, we believe that the issue of order at very 
low energies in the S* = 1/2 KHAFM at (M) = remains 
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a challenging problem. 
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